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^ I We discuss in this paper various aspects of the off-critical 0{n) model in two dimen- 

5^ ! sions. We find the ground-state energy conjectured by Zamolodchikov for the unitary 
minimal models, and extend the result to some non-unitary minimal cases. We apply our 
results to the discussion of scaling functions for polymers on a cylinder. We show, using 
the underlying N=2 supersymmetry, that the scaling function for one non-contractible 
polymer loop around the cylinder is simply related to the solution of the Painleve III 
differential equation. We also find the ground-state energy for a single polymer on the 
cylinder. We check these results by numerically simulating the polymer system. We also 
analyze numerically the flow to the dense polymer phase. We flnd there surprising results, 
with a Cgg function that is not monotonous and seems to have a roaming behavior, getting 
very close to the values 81/70 and 7/10 between its UV and IR values of 1. 
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Since the seminal paper on conformal field theory there has been steady progress 



in understanding two-dimensional critical systems (see 0). Properties at the critical point 
are by now (almost) fully controlled, and new developments have dealt with the vicinity 
of the critical point, where the ultimate goal is to compute all scaling functions. Most 
of the interesting models turn out to be integrable. Factorizabilty then allows one to 
find the ^S-matrix exactly [Q], and use it to calculate the exact ground-state energy g]. 
Although many scaling functions are known in principle, their explicit computation is a 
formidable task This is in contrast with the beautiful and much older case of the 
Ising model, where it was shown P] that spin correlation functions follow from solutions of 
the Painleve III transcendent equation and its generalizations. Recently however, similar 
equations have been met when studying ground-state properties of N=2 supersymmetric 
field theories 0. Using the underlying N=2 supersymmetry in polymers |§] we are able 
here to make contact with these results and to exhibit a scaling function in the polymer 
problem that is given by Painleve III. Remarkably, this property is derived directly from 
continuum considerations, and awaits a lattice derivation. 

We start by considering the 0{n) model. Although results are expected to be universal 
we refer for concreteness to its realization on the hexagonal lattice Although the 
model is defined originally for n integer, the high-temperature expansion allows analytic 
continuation to arbitrary n. The partition function reads 



where the sum is taken over all configurations of self-avoiding and mutually-avoiding loops 
that can be drawn on the lattice, and each elementary edge that belongs to a loop is 
referred to as a monomer. Although this geometrical definition involves global properties, 
locality can be restored if one considers instead oriented loops with a weight exp(+za/6) 
for each left turn, and exp(— za/6) for each right turn. When n = 2 cos a, the relation 
l^left turns — fright turns| = 6 for a closed loop on a planar hexagonal lattice gives the 
correct weight for each loop. (We discuss later what happens on a cylinder.) 

The model has interesting phase transitions for n between -2 and 2. The standard 



(also called dilute) critical point occurs for Kc = [2 -f (2 — n) 2 ] 2 . Parametrizing 
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n = 2 cos — , 
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with 2 < t < oo, the continuum theory at K = K^, has central charge |TC 



c = l-6/t(t+l). 



(3) 



The phase K < Kc is disordered. The phase K > Kc is massless and is referred to as dense 
phase, because the polymers cover a non-vanishing fraction of space even in the continuum 
limit. The continuum limit does not depend on the value of K, and is a conformal field 
theory with central charge 

c= l-6/(t-l)t. (4) 

In the dilute critical theory the thermal operator (i.e. the one coupled to K) is de- 
scribed as $13 with conformal weight h = As discussed in this perturbation 
is integrable. We discuss mainly the case K > Kc in the following. The corresponding 
S- matrix has been studied in |jTl[] ; the integrability requires that the A^-body S-matrix ele- 
ments factorize into a product of two-body ones and thus obey the Yang-Baxter equation. 
One starts by considering the 0(n) model for n integer and assumes that there is a particle 
for each color, whose trajectories may be thought of intuitively as the loops in (J^). The 
0(n) symmetry allows the S-matrix to be written in terms of a set of invariant tensors, 
where these tensors obey an algebra depending on n. Using this algebra, the Yang-Baxter 
equation reduces to a single functional equation with scalar quantities, where n appears 
as a parameter. Analytic continuation becomes possible, although for n non-integer it 
does not define a bonafide S-matrix. However, quantities we calculate using this S-matrix 



should be analytic in n. Therefore, the form-factors calculated in |jTl|] and the ground-state 
energy levels derived here are valid for any n. 



As was also observed in |T^ in a slightly different language, the main property of the 
S-matrix of [l^] is that the invariant tensors obey the Temperley-Lieb algebra [l^] (with 
parameter n). As in lattice models, many interesting properties do not depend on the 
particular representation of this algebra. Thus we can switch from the representation used 
1TT[| , defined only for n integer, to the 6- vertex model representation [jl^, defined for 



m 



any n. In the present context, this amounts to considering the sine-Gordon S-matrix [|T5 
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^Eii + qE22 — Ei2 — -E'21, 
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q = e*^/*, {Eij)ki = SikSji. In this new picture the massive particles are the sine-Gordon 
sohton-antisohton pair whose trajectories can be thought of as the oriented loops in the 
local formulation of the 0{n) model. The S-matrix (Bf) differs from the usual one by a 



gauge transformation that ensures proper algebraic properties [|T6|,|T^ , but our calculation 
is invariant under this transformation. There are no bound states in the spectrum for the 
sine- Gordon couplings we will study. 

We shall derive in the following the ground-state energies E{R) of the 0{n) model on 
a cylinder of radius R and large length L. We will calculate these energies for a variety 
of boundary conditions around the cylinder. It is well known [ p!8| that when one switches 
from the S-matrix to the sine-Gordon S-matrix, boundary conditions must be treated 
very carefully. This is similar to the situation in the lattice model when one switches from 
the geometrical definition to the local one involving oriented loops. On a cylinder, there 
are non-contractible loops in the 0{n) model, and for these the number of left and right 
turns are equal [^. In the local formulation they therefore get a weight 2 instead of the 



desired value n. This boundary condition does not exactly reproduce the 0(n) model; it 
amounts to introducing at infinity an operator of conformal dimension /i = — l/4t(t + l) 

A useful check on the results is to to take the ultraviolet (mass 0) limit of E{R), 
because a general result of conformal field theory [^,22] gives 



Euv^--^{c-12h-12h), (7) 

where h and h are the conformal dimensions of the operator creating this state. The ground 
state is that created by the lowest-dimension operator; only in unitary theories does this 
have h = h = 0. In non-unitary theories one often defines an effective central charge c^q 
as —QREiyR)/^. Thus in the case where non-contractible loops have weight 2, Cgg = 1 in 
the UV. We will show that, as expected, the behavior of E{R) under thermal perturbation 
is given by that for the sine-Gordon model with unmodified boundary conditions, which 
indeed has cuv = 1- 

The technique we will use to calculate the ground-state energies from the S-matrix 
is known as the thermodynamic Bethe ansatz (TBA) [Q. The TEA gives the free energy 
/ = —TlnZ of a particle gas on a line of large length L at a temperature T. This is 
equivalent to considering a spacetime of a cylinder of radius R = 1/T and length L. If 
we change our point-of-view and think of the spatial direction as (Euclidean) time, only 
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the ground state contributes to Z , since only the lowest-energy state can propagate over 
the large distance L. Thus Z = exp(— L£'(i?)), where E{K) is the ground-state (Casimir) 
energy with space a circle of radius K. The equivalence between the two approaches (often 
called modular invariance) means that E{R) = — In Z/L = Rf/L. Introducing chemical 



potentials allows one to change the boundary conditions in the R direction . The TB A 
equations for arbitrary chemical potential ||2^ give the ground-state energy 



^{A}(^) = -Elr / ^^^cosh^ln(l + A„e-^''(^)), (8) 
where the obey the integral equations 

/J at 
—Me-e')Hi + \,e-'^^''^). (9) 
ZTT 



The index a labels the particle species, are fugacities, nia are the particle masses and 
the 4>ab depend on the S-matrix. If the S-matrix is not diagonal, some of the particles in (]|) 
are zero-mass "pseudo-particles" , which result from the diagonalization of the S-matrix. 

To do the TBA calculation, we need to diagonalize the sine-Gordon S-matrix, meaning 
that we find the eigenvalues A(^i|^i, ■ • ■ 5 ^jv) for bringing the ith particle of momentum 



Pi = msinh^i through all the others p5| , |26|] . These eigenvalues are defined precisely 
in the Appendix. Diagonalizing the sine- Gordon transfer matrix is formally equivalent to 
diagonalizing the XXZ model transfer matrix, so doing our thermodynamics is very similar 



to doing the thermodynamics of the XXZ model, which was considered in p^. Putting 



periodic boundary conditions in the L direction quantizes the i-th momentum via 

e^P^^A(^,|^i,02,...,M = l- (10) 

This enables us to define the density of states Po{0) for the particles, whether solitons 
or antisolitons (since the scattering is elastic, the total number of particles is conserved). 
Taking the logarithm of ([1^) gives 

2ml = msinh^i -\- /mln A(^i|^i, . . . , ^^r) 

where n is an integer. Then 

2nPo{e,) = 2n— = mLcoshe, + —Im\nA{e,\eu . . . .On). (H) 



To find the partition function, one takes the thermodynamic trace over aU the eigenstates, 
which is equivalent to summing over aU sohton and antisohton configurations. The diffi- 
culty of this study depends on arithmetic properties of the parameter t. 

We shall first consider the case t integer. In this case, the eigenvalues for the XXZ 
transfer matrix can be written as functionals of densities of pseudoparticles called a-strings 
with a = l,2,...,t— 1 and an antipseudoparticle (1~). Call P+ the associated densities, 
P~ the corresponding densities of holes, and Pa = Pa + Pa ■ Call also Pq the density 
of real particles — the states which are actually occupied by solitons or antisolitons; the 
density of corresponding holes Pq = Pq — Pq . As derived in |2^] in the case t = 2 and in 
the general case in the Appendix, the densities obey the equations 

27^ Pq{6) =mL coshe + (t)i2* pi, (12) 

and 

27rP„ = ^ ^„f, * (13) 
b 

where (j)ab{d) = ^ab(cosh(6'))~^, po = Pq , Pi- = Pi-^ ^"^^ Pa = Pa otherwise; * denotes 
convolution. The sum is over all particles, real and pseudo; lab is the incidence matrix for 
the diagram 

o 1- 

1 t-3 / 

(8) O O O ^-2 

\ t-i 

Each open node represents a string or the antipseudoparticle, while the node labelled ® 
represents the real particle. 

After minimizing the free energy with respect to P^ subject to the constraints ( [T^ ) 
and (|TB|) and defining 

Pa{e) l + A.e-4^)' ^'^^ 

one finds a system of the form (|^) where each node in the above diagram is a particle 
species (the open nodes have = while the node labeled ® is massive), and (f>ab{d) is 
same as above. 
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Taking the UV limit of (^) involves a few tricks [Q. We define Xa = exp(— ea(0)), and 



ya = exp(— ea(oo)). One finds that p3 



CUV = E{mR ^ 0) = — V {Cx, {xa) - Cx, {ya)) 



(15) 



where 



Cx{x) 



'c 



eXe~ 



+ ln(l + Ae-^) 



1 + Ae 

with the contour C going from x to infinity. A change of variables shows that Ci{x) 
L{x/{1 + x)), and C-i{x) = — L(x), where L{x) is the Rogers dilogarithm function 



(16) 



L(x) = ~, 



1 



dy 



Iny ^ ln(l - y) 



y 



It follows from (|^) that the constants the solutions to the equations 



\\{l + \bXb) 



(17) 



where Nab = ItF / d9(t)ab{0). The constants ya in ([T5| ) are nonzero only for those species a' 
with ma' = 0, where they are the solutions to 



ya'=\{{l + \byb')''^'": 



(18) 



where h' also runs only over massless species. These equations for the UV limit are true 
only when (j)ab = 4>ba- 

We have set up the TBA formalism. We only need to specify the boundary conditions 
in the R direction by fixing the fugacities. If we set all Aa = 1, then this corresponds 
to taking periodic boundary conditions: the modified 0{n) model where non-contractible 



loops have a weight of 2. It is easy to find from (p!5|)-(|T^) (and some dilogarithm identities 
|2^ ) one obtains cuv = 1 as expected. 

To get the actual 0{n) model TBA there are two possible approaches. The first 
one is to notice that in the lattice model, non-contractible loops on the cylinder can 



be given a weight n if one introduces a seam |]211]. Using the analogy between particle 



trajectories and the loops, this amounts to computing Tr q-RH ^2i-K-yF ^j^j^ ^ _ iji^ where 
F = i(^solitons — ^antisolitons). More formally, this relation follows from the fact 
that tr e = n for the original 0{n) representation while tr e = 2 in the sine-Gordon 
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replacement, even though both obey the same Temperley-Lieb algebra. From the Bethe 
Ansatz calculation outlined in the appendix, one finds [P7|^ 

F = ll (i^(pi-(^)-p(t-i)W), (19) 
The appropriate trace is therefore reproduced by choosing the fugacities 

Ai- = e-^^^^ A(i_i) = e^^^S (20) 

while A = 1 for the other species. For general 7, this amounts to giving non-contractible 
loops a weight of 2 cos 771. 

In the UV limit, one finds, using (p!5D-(p!8[) (we end up with dilogarithms at complex 
argument) and taking the second derivative of cuv with respect to 7 

For the 0{n) model ground state, 7 = 1/t so we recover (||). Notice that this choice 
corresponds to ^(t-i) = -^1- = —1, while A = 1 for the rest. This ends up removing the 
nodes t — 2, t — 1 and 1~. Indeed, for these nodes, the solutions of (^7[) are the same as 
(PI): Xi^t-i) = Vit-i) = Xi- = Hi- = I for the two end ones, and X{^t-2) = 2/(t-2) = for 
the node joining them. When Xa = Va, ^a{0) = ea(oo), which in fact means that ea{0) does 
not depend on 9. Thus e(t_2)(6') = 00 for all 6*, it no longer contributes to the equation 
for 6(4-3), cind we see that these three nodes decouple from the rest of the system without 
contributing to anything. Thus the 0{n) TBA system at n = 2cos(7r/t) is described by 

12 t-4 t-3 

— o — o o — o 

where all the fugacities A^ = 1. 

Although the 0{n) model at t integer is not identical to the t-th perturbed unitary 
minimal model fl^, their ground states coincide. Thus we have proven that the minimal- 
model TBA is described by the above diagram, as conjectured (and proven for t = 3,4) in 
Another way of deriving the TBA for the unitary minimal models would be to con- 



sider a scattering theory that uses the RSOS S-matrix [30, Iq] instead of the sine-Gordon 



S-matrix |25]. Again, as was discussed in the context of lattice models [18| or integrable 



perturbations [31], the S-matrix obeys another representation of the Temperley-Lieb alge- 
bra, so the thermodynamics deduced from these two cases are equivalent provided one 
considers the sine-Gordon or XXZ system with twisted boundary conditions. The mere 



effect of this is to forbid the 1~ antipseudoparticle and the {t — 1) string [^. The {t — 2) 
string decouples too and we are back to the above discussion. 

The mass m in TBA depends on the 0{n) lattice parameters as m = {^^f^^)''*^^'^^'^: 
where s is a non-universal scale. Notice that for 7 = the same TBA equations would 
describe the ground state energy of say the 8-vertex model. In that case however, if S 
is the Boltzmann weight of the additional vertices that do not conserve U{1) charge, the 
mass m would scale as 5*/^, i.e. with a different exponent. 

We now discuss in more detail the case n = 0, which corresponds to t = 2. As is well 
known this describes the physics of polymers: for small n we the partition function ([^) 
allows only a single loop, the polymer [^. In that case the algebraic properties of the 
S-matrix i in |Tl[] can also be obtained by acting in a Z2 graded space, the loop cancellation 
occurring here due to boson- fermion symmetry. The minimal such S-matrix is equivalent, 
up to another gauge transformation, to the above sine-Gordon S-matrix. 

As discussed in Q the continuum limit of polymers or 0{n = 0) model is described 
by (twisted) A^=2 supersymmetry ||3^. This is true only when n is exactly 0, so only 
properties which can be described precisely at this point are related to supersymmetry. 
The thermal perturbation $13 corresponds in the N = 2 language to the top component 
of the chiral primary field X , hence it preserves the A^=2 supersymmetry. Indeed, 
the S-matrix for this model is precisely the sine-Gordon S-matrix for t = 2 discussed above 



37,26 



The diagram describing TBA in this case has only 3 nodes, one of which is massive. 
With 7 = the ground state energy is minus the logarithm of the largest eigenvalue of 
a polymer transfer matrix that describes non-contractible loops only, with fugacity 2 i. 
The values of Cgg^(mi?) can be extracted from numerical solution of the TBA equations. 
They can also be measured on the lattice by numerical diagonalization of the polymer 
transfer matrix. We considered the honeycomb lattice with cylinders of radius up to 



^ This is true in the regime we are interested in, but not always: see |32] 

^ The S-matrices we refer to here differ from the one in [11 1 by action of the permutation 
operator 

^ Because the lattice bulk free energy vanishes in the polymer case, we can directly identify E 
in the numerical study of c^q and d 
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R = 6 honeycomb faces. Agreement between TBA and lattice computations was found to 
be excellent, even though s must be determined numerically before most comparisons can 
be made. The corresponding curve is given in fig. ^ for R = 6. The main feature of this 
curve is the inflexion point at which the value of c^q extracted from TBA is .7714. Lattice 
estimates of this value are 

R= 2 3 4 5 6 

.732 .748 .754 .758 .762 

which converge nicely to their expected limit. Although c^q is not very interesting as a 
geometrical object, we think this numerical check is useful. In particular the fact that 
thermal perturbation in polymers does not break A^=2 supersymmetry, although reason- 
able, is a little difficult to derive from lattice analysis (due in particular to twisting) and 
is indirectly checked here. 

In the perturbed N=2 theory, a very interesting quantity is Tr e~^^ F{—)^ . It can of 
course be computed from the above by introducing a fugacity 7 close to 1/2 and taking a 
derivative with respect to 7, hence appearing as a solution of an integral equation. What 
is remarkable however is that, as discussed in [^^, it is a kind of supersymmetric index, 
and in this case is related to the solution of a simple differential equation of Painleve III 
type, the radial part of the sinh-Gordon equation. Introduce the function u, the solution 
of 

Uzz + z~^Uz = sinhw (22) 
where u has no singularity on the real axis and has asymptotic behavior as z ^ 

u{z) ~ — 2/31n(^) + constant (23) 

Introduce on the other hand the generating function 

Z^c = 5^i^#"^°''°"^^'' (24) 

where the sum is taken over all configurations of a single non-contractible loop on 
a cylinder of length L and circumference R. One sees from (|l|) and the identification of 
fermion numbers for the 0{n) model and the sine-Gordon theory already discussed that 
this is precisely Tre~^^F(— )^. The results of then give 



Znc{z = mR) ^ —-^^z^u{z) as L/i?—> 00. (25) 
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In particular, at the critical point one has Z^c ~ 3^- The differential equation (p^ ) is 
easy to solve numerically , and for the lattice quantity it can be evaluated using a simple 
modification of the transfer matrix used in the previous problem. Introduce 

d = UuiL/ji^^ ^Znc, (26) 

d goes to 2/7r in the UV limit and goes to zero in the IR. Lattice measurements for R = 5 
are represented in fig. |^. Once again, agreement with (|25| ) is very good. For instance, the 



value of d at the inflexion point is .5104. Lattice estimates are 

R= 2 3 4 5 

.492 .498 .502 .504 

The knowledge of gives access to very non-trivial geometrical information. Call 
11 = 1/ Kc the effective connectivity constant for a given lattice. Then the number of con- 
figurations of a non-contractible loop made of monomers on a cylinder of circumference 
R and length L is expected to scale as 

u;n{R) ~ LR-y^fi^sf{Ns/R'^/^) (27) 

where / is a universal function and s is the lattice-dependent scale factor. From the above 
one finds 

d{mR) = ^ j e^'^^^'^'yf{y)dy (28) 

where m = ( ^sK^ ^^^"^- '^^^ scaling function / is just simply obtained by the Laplace 
transform of the solution of the Painleve equation. 

Another interesting geometrical quantity is the generating function Z defined in a way 
similar to (p^) but summing also over configurations where the loop is contractible. To ob- 
tain the asymptotic behavior of this function is considerably more involved. Indeed, while 
non-contractible loops have a behavior fully determined by the N=2 theory, contractible 
ones depend on properties of the 0{n) model for n close but not strictly equal to 0. We 
therefore have to discuss how to generalize the above analysis to the case where t is close 
to 2. We shall restrict to 

The analysis is easier when t is rational, so we take M to be an integer. This corresponds 
to the non-unitary minimal models (2M + 1, 3M + 1). Relying on the analysis of 
we find that the free energy is as in (P), and the integral equations are of the form (P). 
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There is one massive particle, which we label by as before. There are M + 2 massless 
pseudoparticles. The ones labelled a = 1, M, M +1 are \2a — 1 [-strings (antistrings) if a 



is odd (even), while M + 1 is a 2-string. It turns out to be simpler to write out the TBA 
integral equations (^ in a form where (j)i2 is not symmetric. Then we have 



^^^(e) = j dke 



cosh 6* 

,,,cosh((i-^)f; 



M 

h2{e) = -cP2l{9) 



^ cosh ^ cosh ^ j^2Q^j 



cosh MO 
M 

for j = 2, . . . , M. All the (pab other than (pi2 are symmetric, and all the remaining ones 
are zero. This system has an incidence diagram similar to that at t integer, although here 
particle 1 couples to itself. 

The system with all A = 1 gives of course cuv = 1 for all M, since it corresponds to 
sine- Gordon at different couplings. 

To get the actual ground state energy (which is not the central charge) of the non- 
imitary minimal models, we change the fugacity for the two end nodes, just like in the 
unitary minimal models, by setting Xm+i = ^ m+i ~ while the rest remain 1. This 



again has the effect of merely removing the three end nodes (M, M + 1^ M + 1). We have 
checked that the correct value cuv = 1 — 6/ (2M + 1)(3M + 1) is obtained in that case. 

To observe the actual central charge of the minimal model one needs to project out 
the negative-dimension ground state. To do this, we set also Ai = — 1. As before the three 
end nodes decouple. Let us give some details of the UV calculation. Our choice of (pab 
is not symmetric; symmetric equations are found by rewriting the equations in terms of 
—£2, —£3, ... — cm-i- (We do not change the sign on eo or ei.) The UV free energy is then 
proportional to 



i=2 
11 



The solutions of and (|T|) are 



sin 



xo 



3M+1 



sm 



3M+1 



1 — Xi 



sm 



3M+1 



Sin 



1 + 



sm 



3M+1 

(M+l-j)7r 
3M+1 



sm 



1-2/1 



sm 



2M+1 



sm 



AfTl 



(32) 



sm 



2M+1 

(M+l-j)7r 
2M+1 



3M+1 



sm 



2M+1 



We checked numericaUy that these give the correct answer cuv = 1 — 6M^/ (2M + 1)(3M + 
1). In the M — > oo hmit, the contribution from xq is order 1/M, while the contributions 
of the others are individually However, there are order M such contributions, so 

these also must be included. By replacing the sum over j /M with an integral one checks 
directly that cuv ~ 5/6M as M — > oo and therefore Z ~ 

We can recover the t = 2 case (M = oo) by noticing that except for (poi, all the non- 
zero (pab = TTd{9) at M = OO. Then all the equations (H) except for a = 0, 1 are not integral 



equations any more, but difference equations in M. As detailed in |]23], solving these gives 
the t = 2 result from above. The leading correction involves an infinite number of coupled 
integral equations, which we have not been able to simplify noticeably. Numerical solution 
for finite M and extrapolation to M infinite produces results in good agreement with the 
lattice computation. Results in the latter case for 



ttR 



(33) 



are represented in fig. |^ for a cylinder of R=5 honeycomb facesS. 

The flow from the critical to the dense 0{n) model should also present interesting 
features, but the scattering theory is not so well-understood and we postpone most of its 
study to a subsequent publication. We shall content ourselves with the following observa- 
tion. For unitary theories, one can define a c-function, which is the central charge at the 
conformal point and which never increases under renormalization group flows (increasing 
mR) . It has been conjectured several times that the ground-state energy on a cylinder 
plays the role of an alternate c-function, even for non-unitary theories. The dilute- and 
dense-polymer theories are non-unitary. Moreover, in the sector where non-contractible 



In this case a non-universal bulk free energy term must be subtracted to get E. 
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loops get a weight 2, they both have c^q = 1 in the UV. If a Cgg-theorem held, Cgg(m-R) 
here should be a constant in the flow from dilute to dense. We have estimated this quantity 
numerically using our polymer simulation. As always, in the dense region the convergence 
is a little worse than in the dilute region. But the main features of the resulting curve 
proved very stable, with uncertainties of only around 1% for the asymptotic quantities. 
The curve for R = 6 hexagonal faces is given in flg. |2| . It violates the c^q = 1 theorem in 
a spectacular fashion. Moreover, the maximum and minimum lie respectively very close 
to 81/70 and 7/10. These are the two values of the central charge neighboring c = 1 in 
the minimal = 1 series. Recall that the A^=2 theory related to polymers belongs also 
to the A^=l series. The flow between dilute and dense polymers seems therefore to be an 
example of a "roaming trajectory" as introduced in HO . 
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Appendix A. Derivation of sine-Gordon TBA equations 



In this appendix, we derive the equations (^2]) and (|T^) for the density of states of the 
sine-Gordon model at integer t. The allowed A{9i\9i, 62 - ■ ■ 6* at) in ([Tl| ) are the eigenvalues 
of the sine-Gordon "transfer matrix" Tah{0i\6i, 62 ■ ■ ■ On) for bringing the i-th particle (with 
rapidity 9i) through the other — 1 particles and ending up with itself. The components 
of Tab{Oi) can be written in terms of the S-matrix elements as 

iTat{e.))\'!^l = st%io - ^iXt -02)... si:'_:^je - e^l (A.i) 

{fc,}=0,l 

where the cj (dj) label whether each of the A^ particles before (after) the scattering process 
is a soliton or antisoliton; the kj are the intermediate particles. The sine-Gordon S-matrix 
is given in (^), where 
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Formally the sine- Gordon "transfer matrix" is equivalent to the transfer matrix of 
the six-vertex model or of the XXZ spin chain. To diagonalize it, we can rely on the 
well-discussed technique of the Bethe ansatz The only additional ingredient here is 
the prefactor Z{9); by the definition (|A.1|) we see that its presence causes the eigenvalues 
A{9i) to be multiplied by Y[^=i ^i^i ^ ^j)- The result for integer t (for example, see 
or p^, or [|6l for the case t — 2), is that 

^ lnA(^,) = ( A/^inZ) * P+ - * P+ * P+ , (A.3) 



and 



(A.4) 



Pi- = Pit-i) 

where the sums run over the strings (but not the antipseudoparticle) . The kernels in these 
equations are defined most easily in terms of their Fourier transforms: 

B (k)- f ^^^^^B (0) - S o ^o^hf sinh(t-a)f sinh6f 
0„,(/c)-y-e 0^,(6) -5., -2 ^j^^f sinhtf ^^'^^ 



for a > 6 with Ba& = B^a, and 



_ sinh(t-a)f 
^"^^^ " sinhtf 



(A.6) 



Pi- ^P{t-i) 
Qb = 0(t-i),b 

To simplify these equations, it is much more convenient to work in Fourier space. We 

d 

d8i 



define z{k) as the Fourier transform of -A-Im\rvZ , so that 



Tik sinh(t-l)^ 

2z{k) cosh — = ; — — = pi 

^ ' 2 sinh^t ^ 

It is easy to prove the following trigonometric identities: 

2Qab cosh — = Qa+l,b + Qa-l,b " Sa+l,b - Sa-l,b 

ttIv 

2pa cosh — = Pa-1 + pa+l (A-T) 
2qa cosh — = qa-l + Qa+l - Sa,t-2 " Sa,f 
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Using these in ( |A.4|) and ( |A.3| ) gives 

2Pa cosh ^ - Pa+1 - Pa-1 = - + Sa,t-2P+ 

2P(,_i) cosh ^ - P(,_2) = -P+_2) (A-8) 

2[(/c) cosh — - Pi = -P+. 
2 

where Z(/c) is the Fourier transform of ^/mlnA(^), and a = — 2. Using P~ = 

P — P+ and the fact that the Fourier transform of 1/(2 cosh ^) is 1/ cosh^ gives us ( |T^ ) 
and (|T3|). 

From the Bethe ansatz one also finds that the number of sohtons minus the number 
of antisoUtons is 

I d6(^-P+{e) + 2Y,aP+ie) + 2P+{6)^ (A.9) 

In the XXZ language, this is twice the value of the third component of the spin. Using the 
relation ( |A.4| ) for a = p — 1 gives the relation (p!9|). 
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Figure Captions 

Fig. 1. Lattice measurements for the various scaled ground states c^q: non-contract ible 
polymers only, each with a weight 2 . d: one non contractible polymer, dd (c): 
one polymer, contractible or not. 

Fig. 2. Flow of CqQ to the dense phase. 
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